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lar viscous fluids driven by a genuine cylindrical fractional Bronwnian motion with 
Hurst parameter H € (4,2) un der Dirichlet boundary condition on 2D square domain. 
First we prove the existence and regularity of the stochastic convolution corresponding 
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dynamical system generated by non-Newtonian fluids has a random attractor. 
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1 Introduction 

In this paper, the stochastic non-Newtonian fluids driven by fractional Bronwian 
motion (fBm, for short) on [0, n] x [0, n] are studied. The constitutive relations for 
such fluids were introduced by Bellout, Bloom and Necas [1] to describe the isothermal, 
nonlinear, incompressible bipolar viscous fluid. It has the form 

T ij = -pSij + 2/x (e + | e | 2 ) ae^ - 2/xi Ae^ (1.1) 

n Jk = 2^ k (1.2) 

where t%j is the components of the stress tensor, t^ is the components of the first 
multipolar stress tensor, and p is the pressure, eij are the components of the rate of 
deformation tensor, i.e. 

««=K£ + £)- (L3) 

e,fio,fii > 1 and a, < a < 1, are constitutive parameters. The constitutive relation 
(jl.ip and (jl.2p . and the condition of incompressibility, yield the following nonlinear 
partial differential equations (we call it Bellout-Bloom-Necas fluids): 

du 
Plu + P( u • V)u + Vp = V • <jt(u)e - 2//i Ae) + pf (1.4) 

V-« = 0, (1.5) 

where p is the constant density, n{u) = 2/xo(e + |ej 2 )~ "2 is a nonlinear viscosity, and / 
is the external body force vector. 

There are many works concerning existence and regularity of solution to the Bellout- 
Bloom-Necas fluids and its dynamics (see, for instance, [2j [IJ [31 HI [5j |6] ) . In this paper, 
we consider the stochastic Bellout-Bloom-Necas fluids driven by a genuine cylindrical 
fBm with Hurst parameter H G (4,5): 

du , „. „ „ , , , „ . , dB H (t) 
— + (u ■ V )u + Vp = V • (ji(u)e - 2/iiAe + / ; 
at at 

u = 0, T ijk r]jrj k = 0, 
u = u , 



x£(D, t > 


(1.6) 


ieO, t > 


(1.7) 


x G SO, t > 


(1.8) 


iG0,i = fl 


(1.9) 



where O is a 2D square, i.e., O = {[x\,X2} | < X\ < ir, < %2 < vr}. The fractional 
Brownian noise enters linearly in the equation and the fBm models the noise source. 
The fBm is a family of Gaussian processes and some useful properties of these process 
were given by Mandelbrot and Van Ness [7] . For H < ^ the fBm is not a semimartingale 
and they can be used in modeling phenomena with intermittency and anti-persistence 
such as financial turbulence. 

The preprint [8J treats stochastic Bellout-Bloom-Necas fluids that the noise term 
has a trace-class correlation, and moreover they treat the case H > ^, which allows 
one to solve the equation using stochastic integrals understood in a pathwise way. In 
this paper we first provide a detailed study of the existence and regularity properties 
of the stochastic convolution corresponding to stochastic non-Newtonian fluids driven 
by fBm with H € (■?, ^)- The approach for dropping the Hilbert-Schmidt operator 
hypothesis follows [9 J . Then we establish the existence of solution by a modified version 
of fixed point theorem in some specific Banach space. The a priori estimate for solution 
in intersection space follows [HJ. Finally we construct the random dynamical system 
associated by non-Newtonian fluids and prove the existence of random attractor follows 
the framework of |10j . 

We emphasize four points in our paper, (i) By careful estimation, the growth 
speed for eigenvalues of the linear differential operator ensures the convergence of the 
stochastic convolution respect to fBm integral. When investigating the dynamics of 
Bellout-Bloom-Necas equation perturbed by fraction Brownian noise, (ii) the ergodic 
property of stochastic integral with respect to fBm ensures us to construct a random 
dynamical system, (iii) the at most polynomial growth of sample path ensures the 
existence of random absorbing sets, and the regularity of stochastic convolution with 
respect to infinite dimensional fBm with H E (j, |) ensures the existence of absorbing 
set with higher regularity, (iv) Under certain condition ( which limits the injection 
parameters of fractional integral space to 1? space ) , the random dynamical system has 
a random attractor as the case H > i in [8j. 

The rest of the paper is organized as follow. In section 2, we formulate the mathe- 
matical setting for stochastic Bellout-Bloom-Necas fluids with Dirichlet boundary con- 



dition on 2D square domain and recall the Wiener integrals with respect to infinite 
dimensional fBm by the framework of [9]. Section 3 is devoted to the existence and 
regularity properties of the stochastic convolution. In section 4, the global existence 
and uniqueness of solution is obtained. In section 5, we prove the existence of a random 
attractor for the random dynamical system generated by non-Newtonian fluids. 

2 Preliminaries 

We use the standard mathematical framework of this model. 
V = U= (4> x , <j) 2 ) € (C °° {O)) 2 : V • <p = and <f> = on do\. 
H = the closure of V in [I? {0)j with norm • |. 



V = the closure of V in (ff 2 (C)) with norm 



v ■ 



H k = the closure of V in (H k (0)) with norm | ■ \j. for k £ N. 

For other Banach space Y, denote the norm | • |y. 

B Y (M) :={\x\ Y <M,x£Y}. 

Let H' and V' be the dual spaces of H and V respectively. It follows that V C H = 
H' C V and the injections are continuous. 

Denote by (■, ■) the inner product in H, < -, • > the dual pair between V' and V. 

C(H,V) := { all the bounded linear operator from H to V }, C{H) := C(H,H). 

C2{H) := { all the Hilbert-Schmidt operator from H to H }. 

We use lowercase Ci,i E N for global constants and capital C for local constants 
which may change value from line to line. 

Inspired by [II], we select space X := C([0,T];H) f| L 2 (0,T; V) with norm \-\x-= 
I ' \c([o,T];H) + I ' \l 2 (o,T;V) f° r solutions. For the completeness of space X we refer to [8]. 

First we define a bilinear form a(-, ■) : V x V — >M., 

a(u,v) = -(Au,Av). (2.1) 

Proposition 2.1. The rate of deformation tensor e(u) and bilinear form a(-, •) have 
the following properties: 



(i) 

V • e(n) = -An Vn E V, (2.2) 

V • (Ae(n)) = ^A 2 n Vn E if 4 . (2.3) 

(mJ For u, v E V, 

1 2 £ fde^de^ * fde^de^ul 

2 V ; .^7 cb,- 9x fc .^ / &E fc 5x fc v ; 

J,j,fc=l i,j,k=l 

(Hi) (f^ Lemma 2.3) There exist ci,C2 > 0, s.£. 

ci|n|y < a(u,u) < C2\u\ v , Vn E V. (2-5) 



Next we define the abstract differential operator A and the analytic semigroup 
S(-). According to (hi) of Proposition 12.11 we can use Lax-milgram Theorem to define 
AeC(V,V): 

< Au, v >= a(u, v) V«,dEV. (2.6) 

And we have 

Proposition 2.2. \8j 

(i) Operator A is an isometric form V to V . Furthermore, let D(A) = {u E V : 
a(u,v) = (f,v),f E H}. Then A E C(D(A),H) is an isometric form D(A) to H. 

(ii) Operator A is self-adjoint positive with compact inverse. By Hilbert Theorem, 
there exist eigenvectors {ej}^ 1 C D{A) and eigenvalues {Xi} ( *L 1 s.t. 

Aei = X ei , e t eD(A), i = l,2,--- (2.7) 

< Ai < A 2 < • • • < Xi < ■ ■ ■ , lim Xi = oo. (2.8) 

i— >oo 

And {ei} ( £L 1 form an orthonormal basis for H. 

(Hi) For u£ D(A), 

Au = V-(Ae(n)) = -A 2 n, (2.9) 

i.e., A = PA 2 where P is the Leray projection operator from L 2 (0) to H. 



Noticing that A is a self-adjoint positive linear operator with discrete spectrum, we 
define the fractional power of A by following the framework of Chueshov [12] section 
2.1: 

Definition 2.3. For a > 0, 

( oo oo ^ 

D(A a ) := \ h = Y,c k e k G H :J2 c l^kf < oo L (2.10) 



fc=l fc=l 



{oo 
formal serial > c^e^ suc/i £/jo£ > c fe (A^) < oo 
fc=i 



(2.11) 



A Q /i = £c fc A£e fc , kD(n (2.12) 

fc=i 

Let S^q, = D{A a ). Then 3~ Q is a separable Hilbert space with the inner product 
(u,v)j a = (A a u,A a v) and the norm ||w||j a = |j4 q- u|; 3~_ a denotes the union of bounded 
linear functional on 3~ Q ; Particularly, we have 3~o = H, $\/2 = V, 3~-i/2 = V and the 
norm of V and H 2 (0) are equivalent; For a\ > o~i, the space 3" c - 1 is compactly embedded 
into 3~(j 2 . 

Since j4 is a densely defined self-adjoint bounded-below operator in Hilbert space H, 
we deduce that A is a sector operator and it generates an analytic semigroup S G £{H) 
(see, for instance, [13] section 1.3), 

/"OO 

S(t) := e~ tA = / e~ tx dE x . (2.13) 

For a survey of the properties of the analytic semigroup we refer to [8]. 

In order to write down the abstract evolution equation, we need to handle the 
nonlinear terms. Following the method in dealing with Navier-Stokes equation, we 
define the trilinear form: 

b(u,v,w) = Y] / u i7 p- Wj dx Vu,v,w € H%(0). (2.14) 

iJ=1 Jo JXi 

Since V C Hq(0) is a closed subspace, &(■, •, ■) is continuous in V x V x V. From [2] 
we have, 

b(u,v,w) = —b(u,w,v), b(u,v,v)=0 Vu,v,w,£ H (O). (2.15) 

For u,v &V, define the functional B(u,v) € V: 

< B(u, v),w >= b(u, v, w) Vw G V (2.16) 



and denote B(u) := B(u,u) G V. 
For u G V, define N(u) as 



< N(u),v >= n(u)eij(u)eij(v)dx Vt; G V. (2.17) 

Then iV(-) is a continuous from V to V and 

< iV(u), v>=- (V • (/*(u)e(u))) • ucte. (2.18) 

Jo 

Comprehensively, we have the following abstract evolution equation from problem 

(ILED-OD: 

du + (2fi 1 Au + B(u) + N(u)) dt = dB H (t), 

(2.19) 

u(0) = u . 

Since the derivative of fBm exists almost nowhere, we will give a mathematical interpret 
of above equation in the next subsection. Without loss of generality, we set \i\ = 1 in 
the sequel. 

Next we introduce the Wiener-type stochastic integral with respect to fBm. For all 
T > 0, let j3 (t) be the one-dimensional fBm. In this paper we only consider the case 
H < g following the [9] and for a survey of Winer-type stochastic integral we refer to 
|15j . By definition j3 is a centered Gaussian process with covariance 

R(t, s) = E(p H {t)p H {t)) = -(t 2H + s 2H -\t- s\ 2H ). (2.20) 

And j3 H has the following Wiener integral representation: 

f3 H (t)= I K H (t,s)dW(s), (2.21) 

Jo 

where W is a Wiener process, and K H (t, s) is the kernel given by 

K H (t, s) = c H (~\ 2 (t - s) H ^ + s^Fi 1 -). (2.22) 

ch is a constant and 

F(z)=c h (--h) f Z r H -l (l-(l + r) H -^dr. (2.23) 

By ([2321 we obtain 

^<t, 8 ) = c H (H-l){t-s)*-i(?) h - H . (2.24) 



Denote by £# the linear space of step function of the form 

n 

¥>(*) = J>l (ti>tl+l] (t) (2-25) 

i=l 

where n € N, a, € R and by IK the closure of £# with respect to the scalar product 

<l[o, 1 ],l[M>5(=«(M) (2.26) 

For (/? 6 £# we define its Weiner integral with respect to the fBm as 

[ T cp(sW H (s) =f>(&? +1 -Pg). (2.27) 

The mapping 

n . T 

<P = J2 a * 1 (U,t»i]-+ V>(s)d/3 H (s) (2.28) 

is an isometry between &h and the linear space span{/3 (£),0 < i < T} viewed as 
a subspace of L 2 (0,T) and it can be extended to an isometry between "K and the 
span L {f3 H (t),0 <t< T}. The image on an element ^ € IK by this isometry is called 
the Wiener integral of \l/ with respect to /3 . 
For every s < t, consider the operator K* 

(K*)<p(a) = Kit, sMs) + J\<p(r) - ^))^(r, s)dr. (2.29) 

We refer to [16] for the proof of the fact that K* is a isometry between 'K and L 2 (0, T). 
For H < 2 , the reproducing kernel Hilbert space % can be represented by the fractional 
integral space. Namely, 

U = {K* H )'\L\Q,T)) = 4: H (L 2 (0,T)), (2.30) 

where II is the family of functions / that can be represented as a fractional 

IZ- -integral of some function <\> £ L 2 (0,T). As a consequence, we have the following 
relationship between the Wiener integral with respect to fBm and the Wiener integral 
with respect to the Wiener process: 

ip(s)d(3 H ( S ) = f {K* t ip){s)dW{s) (2.31) 

o Jo 

for every t < T and (p € "K if and only if K%<p € L 2 (0,T). Since we work only with 

Wiener integral over Hilbert space, we have that if u £ L 2 (0,T;H) is a deterministic 



function, then the relation (|2.31|) holds and the Wiener integral on the righthand side 
being well defined in L 2 (ft; H) if K*u belongs to L 2 (0, T; H). 

In the following we concern the infinite dimensional ffim and stochastic integration 
(see, e.g. |17j). A standard cylindrical fractional Brownian motion is defined now. 

Definition 2.4. fiffl Let (£1, J 7 , P) be a complete probability space. A cylindrical process 
< B , • >: Q, x R + x H — v R on (Q, J 7 , P) is called a standard cylindrical fractional 
Brownian motion with the Hurst parameter H £ (0, 1) if 

1. for each x £ iJ\{0}, ynqr < B H (-),x > is a standard scalar fBm with Hurst 
parameter H; 

2. for a, j3 £ M. and x,y £ H , 

<B H (t),ax + Py>=a<B H (t),x>+(3<B H (t),y> P-a.s. (2.32) 

For H = g, this definition is the usual one for a standard cylindrical Wiener process 
in H. For the complete orthonormal basis {e n } nG N of H (which is generated by linear 
differential operator A), letting /3„ (t) =< B H (t), e n > for n £ N, the sequence of scalar 
processes {/3^} ng N is independent and B can be represented by the formal series 

oo 

B H (t) = Y J ^(t)e n (2.33) 

n=\ 

that does not converge a.s. in H. Although for any fixed t the series f|2.33j) is not 
convergent in L 2 (0 x H), we can always consider a Hilbert space U\ such that H C U\ 
such that this inclusion is a Hilbert-Schmidt operator. In this way B given by f|2.33j) 
is a well-defined U\ -valued Gaussian stochastic process. 

Let <&(s),0 < s < T be a deterministic function with values in C2(H), the space of 
Hilbert-Schmidt operators on H. The stochastic integral of $ with respect to B is 

defined by 

ft °° ft °° ft 

/ <S>(s)dB H (s) = V / $(s)e n dP*(s) = V / (K*($en))(s)dp n (s) (2.34) 

J ° n=l J ° n=l J ° 

where f3 n is the standard Brownian motion. However, as we are about to see, the 
stochastic linear additive equation in its mild form can have a solution even if L &(s)dB (s) 
is not properly defined as a if -valued process. 

9 



3 Linear stochastic evolution equations and stochas- 
tic convolution with fBm 

In this section, we will work with a cylindrical fBm B on the real separable Hilbert 
space H . First we prove the existence of solution for linear stochastic evolution equation 
driven by fBm with H G {\,\)- Consider the equation 

dZ = AZdt + dB H , Z(0) = u G H (3.1) 

As noted in [S] and |18| . the stochastic integral L IddB H (s) is not well-defined as a 

H -valued random variable since the identity operator Id £ C2{H). We then consider 

the mild form of the equation, whose unique solution, if it exists, can be written in the 

evolution form 

Z(t) = S(t)u + f S(t- s)dB H (s). (3.2) 

Jo 

Remark 3.1. In [18] they consider the noise term as QdB H and assume that <3? G £2 
or S(t)& G Ci- In such case the infinite dimensional fBm noise and be viewed as finite 
dimensional since the Hilbert-Schmidt operator is compact. By contrast the noise we 
consider (i.e., $ = Id) is rougher and there is no reason to assume that $ G £2- 

In order to obtain the existence of the stochastic convolution above, we need the 
follow estimate about the spectrum of operator A. 

Lemma 3.2. The eigenvalues of operator A satisfy 

A m „ > (m 2 + n 2 ) 2 , m,nGN. (3.3) 

Proof It is well known that, the eigenfunctions and eigenvalues of Dirichlet-Lapacian 
on a 2D square (e.g. see [IS]) are 

(f)(xi,X2) = C sin mx\ sin nx2 

(3.4) 

Imn = rn + n , with m, n G N. 
Rewrite the the index of 7 and we get the spectrum 

<*D = {-yj}?=i, 0< 7 i<72<--- (3.5) 

10 



Let as be the spectrum of the stokes operator (i.e. —PA where P is the Leray 
projector) with homogenous Dirichlet boundary conditions (which we refer to classical 
Stokes operator), and we have 

*S = $j}jLi, 0<A!<A 2 <--- (3.6) 

According to Theorem 1.1 in [20], we have 7^ < A& for all positive integers. 
Back to the operator A of non-Newtonian fluids, we have 



PA 2 e, 



-PA(-PA) ej 



X fa 



For technical reason we need the lemma: 

Lemma 3.3. The following inequality holds for all A > 0. 
-A / r x „ \ 2 



Proof. Calculate the integration and we have 

i-X 
e- 2x 



j X ( e v - l)y H -'idy] dx 
e~ 2x I / {e y - l)y H -*dy + (e y - l)y H ^dy dx 



(3.7) 



Thus we can reorder the index and estimate the spectrum of A 

A mn = A^ n > 7^ n = (m 2 + n 2 ) 2 . (3.8) 



□ 



/ e ~ 2X ( F^ V ~ l )y H ~^ d y) dx = C(H)< 00. (3.9) 



A / ,1 , ,max{l,.} N 2 ^ A ^ 



11 



Firstly we estimate the first term in the integral by parts of integration. 



\ e y-l)y H -ldy 



^t(V -l)y H -s[- j\vy H -*dy 



H - 


i 

2 


1 




H - 


1 

2 


1 




H - 


1 
2 



(e — 1) — lim 



e y -1 



o v h- H H-h H + h 



e y dy H+ 2 



2 -/0 



'e-D 



1 



#- J # + 2 



L_ f^* x 



e y y H+ ^dy 



=C{H) < oo. 
Secondly we estimate the other term 



3 



A / ™ax{l,i} 

e" 2x / (e y - l)y H ~'idy \ dx 



A 



/ (e»-l)-i/ H- idy) da; 

( e » _ i) 2 . y 2H ^ e - 2x (x - l)dydx 
(e y - l) 2 • y 2H - 3 e- 2x {x - l)dxdy 



/A rx 
I ( e V - . n2., ( 2H-3„-2x, 

< /" f (,y - d- . „2"-\.-2..-, 

/A /*cx) 

(e y - l) 2 • 2/ 2if " 3 / e" 2 *^ - l)da;dy 



l -(y-l)e~ 2 y + \e- 2 y^dy 



<C I (1-e 



-^2 . y 2ff-2 + (1 _ e - w) 2 . y 2ff-3 dy 



<C/ x ^"W^ 



=c 



1 1 

+ 



2# - 1 2H - 2 



C(fl") < oo. 



Finally we have 



<, < 



-2x 



{e y_ l)y H- 2dy + 

Jl 



max{l,x} 



/•A / /•max{l,cc} . 

< / e~ 2x C(H) + (e y - l)y H ~^dy dx 



(e y -l)y H -2dy\ dx 

\ 2 



A 



A 



<2 / e- 2x C(H) 2 dx + 2 / e" 2:c 



/•oo 

<C(H) / e~ 2a; cix + C(fl") < oo. 
Jo 



max{l,i} 



(e y -l)y H -%dy) dx 



(3.11) 



(3.12) 



(3.13) 



12 



□ 



Remark 3.4. There is a mistake in the proof of Lemma 2 [9j. That is (in Page203), the 
quantity 

POD / px \ 2 

K A = J U (ey-l)y A - l dy\ dx (3.14) 

is infinite. Therefore the theorems which use this lemma (such as Theorem 1.1 in [9] 
and Theorem 2.1 in |21j ) may not hold. However, the argument there can be modified 
by using our lemma and the expected result can be obtained. 

The statement of existence theorem for stochastic convolution follows. 



Theorem 3.5. If the Hurst parameter H > j, the stochastic convolution J Q S(t 
s)dB H (s) is well defined. 



Proof. The proof is base on Theorem 1 of [9]. It is sufficient to estimate the mean 
square of the Wiener integral of (|3.2|) . 



E 



ft 2 °° ft 

/ S(t - s)dB H (s) = E V / S(t- s)e n dp*{s) 

JO H n=1 J0 



(3.15) 



H 



Using (12.29h and the representation (12.31J) . we have 

E\z(t)-S(t)uo\% 
<2^y \S(t-s)e n \ 2 K 2 {t,s)ds + 2Y J f J (S(t - r)e n - S(t - s)e n )^(r,s)dr 



'~-h + h 



By [22] Th3.2, we have 



ds 



(3.16) 



K(t,s)<C(H)(t-s) H -2 S 1 



(3.17) 



Then, 



h<C(H)Y] f\e-^- s)X -e n \ 2 (t-s) 
Z Jo 



2H-l s 2H-l ds 



C{H)Y^{2X n y 2H [ 2Xni 'e-^ H -\t-^-f H -Hv 
„ Jo 



2Ar 



(3.18) 



13 



where C(H) depends only on H and we use the change of variable t — s = -^—. Since 

f2X n t 



e ~ v v 2H - 1 



"-i' 2 "- 1 



dv 



< 



X n t 



t 



o 



e -v v 2H-l { ^2H-l dv+ e -v 2Xnt 

— J X n t 

Xnt 



2A T , t 



2H-1 



(t 



\2H-1 



L\ r . 



dv 



<C(t,H) I < r 

'o 



v„,2H-\m, i (n\ +\2H- / ,-(2A„i-V) ' V 



dv + (2A„t) 
<C(t,H)r(2H)+t" 

<C(t,H) + C(t,H)e- Xnt I 

Jo 



2X n 



2H - l dv' 



■211-1 l at e -K n t e v> v ,2H-l dv , 



(3.19) 



rXnt 



v nH - x dv' 



<C(t,H)-e- Xnt (\ n t) 2H 

<C(t,H), 

we have 



h<C(H)-J2\- 2H -C(t,H) 



<C(t,H). £ 



,-2Uff 



(3.20) 



faP + PV 

< C(t,H) ■ 2(3 D {AH) • ((AH) < oo 
where /3d (s) is the Dirichlet beta function and £(s) is the Riemann zeta function (for 
definition see 



For the second sum, we have 



dK 



l2 = 2 Y,[ f ( 5 (* - r ) e « - S ^ - s ) e «)^( r •' )r/r 



ds 



H 



< 



^E 



f e -(*-r)A„ _ e 






(3.21) 



ds 



By the fact that -g=-(r, s) < for every r, s G [0, T] and 



<9r 



r,s. 



< C{H) -{r-s 



,H- 



(3.22) 



we have 



h < C(H) J2 [ ( [ (e~ (t ~ r)An - e- (t - s)A ")(r - s) H -%dr\ 
<C(H)^2 [ ( [ S ( e -^- s - v)x " - e- {t - s)x ")v H -'Uv] ds 



ds 



WO 



(3.23) 



< 



C(H)Y^ I ( [ U (e-( u - v)Xn 



-u\ n \„,H 



)v zdv) du 
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By the change of variables v = 4f-, u = v - , we have 

A n A n 

< C{H) VJ X- 2H £ nt e~ 2x QfV - l)y*-*dy) 2 dx 

< C7(fl" ) ^ \- 2H C{H) (by lemma £3D 

n 

<C(H)-2p D (AH)-((4H) <oo as 4ff > 1. 



(3.24) 



D 



In the following we investigate the regularity of solution. First we show that B H € 

c([o,oo) ; y')- 

Since 

oo oo 

n=l n=l (3.25) 

= 2/3^(2) • C(2) < oo, 

A~z : H — >• V' is a Hilbert-Schmidt operator from i? to V. Thus the genuine cylindrical 
fBm B can be viewed as a F'-valued fBm with incremental covariance operator A~2 . 
And we have that B has continuous trajectories on V'. 
Let 

z(t) : = I S(t-s)dB H (s), (3.26) 

Jo 

Y(t) :=[ S(t- s)B H {s)ds, (3.27) 

Jo 

and we have 

Lemma 3.6. For H > \, Y(-) belongs to C 1 ([0, oo); V) P-a.s. and 

= B H {t) + A I S(t- s)B H (s)ds ( 3 - 28 ) 

Jo 

= *(*)■ 
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Since the proof is a modification (consider the base space as V) of Lemma 5.13 in 
and Proposition 3.1 in |25j . we only sketch the main ideas. 
Step I: 
When Aq : V — > V is a bounded linear operator, we have 

t 



z(t)= / A z(s)ds + B H (t), t>0,P-a.s. 
Jo 



(3.29) 



Step II: 

Let A n = n(nl — A)^ 1 A, n 6 N (n large enough) denote the Yosida approximations 
of A. Denote by S n {-) the semigroup generated by A n in V' and set z n {t) = j Q S n (t — 
s)dB H (s). By using the so-called factorization method, we have for q > 2 

E sup \z n (t) — z(t)\y, — > 0, as n — > oo. (3.30) 

te[o,T] 

Step III: 

Since z n {-) is a strong solution to the stochastic differential equation in V' 

dz n {t) = A n z n (t)dt + dB H (t), z n (0) = 0. (3.31) 

We have 

z n {t)= fo t A n z n (s)ds + B H (t). (3.32) 

Setting Y n (t) = J Q z n (s)ds, we have Y n {-) is the solution to the initial value problem 

-Y n (t) = A n Y n (t) + B H (t), Y n (0) = 0. (3.33) 



10 

d 



dt 



Thus we have 



and so P-a.s. 



Y n (t) = I S n (t - s)B H (s)ds, (3.34) 

Jo 



z N (t) = A n Y n (t) + B H (t). (3.35) 



clearly 



lim Y n (t) = Y(t) = [ S(t- s)B H (s)ds. (3.36) 

n->oo J 

By step I and step II, we have 

lim A n Y n (t) = lim Ail - -A^Y^t) = z(t) - B H (t). (3.37) 

n— too n^roo n 
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Since the operator A is closed, we conclude that Y(t) € V and AY(t) = z(t) — B H (t) 
P-a.s.. 

The following proposition is an immediate consequence of above lemma. 

Proposition 3.7. If H £ (i, i). Then for alluo E V the process Z(t,x) = S(t)uo+z(t) 
has an V- continuous modification. Moreover, 

z(t) = A[s(t-s)B H (s)ds + B H (t), t>0 (3.38) 

Jo 

holds with probability one. 

4 Solution of stochastic non-Newtonian fluids 

We interpret equation (J2.19P as an integral equation and define the solution as 
follow. 

Definition 4.1. The solution to equation (J2.19P is defined as a function u € C([0, T]; H)f\ 
L 2 (0,T;V) s.t. the following integral equation holds P-a.s.. 

u{t) = S(t)u - I S(t-s)B(u(s))ds- I S{t-s)N(u(s))ds+ f S(t-s)dB H (s) 
Jo Jo Jo 

(4.1) 

where the first and second integral are Bochner integral, the last integral is a stochastic 

integral defined in previous section. 

We seek the solution of 2D stochastic Bellout-Bloom-Necas fluids by the fixed point 
in space X = C([0,T];H) n L 2 (0,T;V). Firstly we will prove the local existence and 
uniqueness results. 

For ii£l, let 



Ji(«) : = - f S(- - s)B(u(s))ds, (4.2) 

Jo 



j 2 [u) ■ = - S(-- s)N(u(s))ds. (4.3) 

Jo 



We have following estimates. 
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Lemma 4.2. [8] J\ : X — > X and for all u, v 6 X, we have 

\Ji(u)\ 2 x < Cl \u\%, (4.4) 

I 7 / \ t / \|2 <s ( \ |2 I |2 | I |2 I |2 \ 2 

KlW ~ J lWIX S C 2 \\U\c([0,T\;H) " I U Il 2 (0,T);V + \ v \C([0,T];H) ' l U lL 2 (0,r);Vj 



(4.5) 



!«-«&• 



Lemma 4.3. ( L 8] J 2 '■ X — > X and for all u, v G X, we have 

\Mu)\x < <%\ u \h(o,T;V)> ( 4 - 6 ) 

\Mu)- J 2 (v)\ 2 x <c A T\u-v\ 2 x . (4.7) 

We apply the following version of the contraction mapping theorem. 

Lemma 4.4. (\2J$ Lemma 15.2.6) Let F be a transformation from a Banach space E 
into E, <p G E and M > a positive number. If F(0) = 0, \</>\e < 2 M and 

\F{u)-F(v)\e<^\u-v\ e Vu,v£B e (M), (4.8) 

then the equation 

u = (j) + F(u) (4.9) 

has a unique solution u G E satisfying u G Be(M). 

We now prove the main result of this paper. 

Theorem 4.5. (local existence and uniqueness of solution) For all uq G H , there 
exists To > s.t. equation (J2.19J) P-a.s. admits a unique solution u G C([0, To];i7) PI 
L 2 (0,T ;V) in the sense of (J4TTJ) . 

Proof. Fix uj G Q. Let 

<P(t) = S(t)u + z(t). (4.10) 

By the properties of semigroup S(-) and Lemma [3.7[ S(-)uo,z G C([0,T]; V) C X. 
Then we have 

\<P\x < \S(-)u \ x + \z\ x < 2|u | + \z\ x . (4.11) 

Let M{u) = 2(2|n | + \z(w)\x). 



II 



Construct the mapping T = J\ + J2, then for all u, v € X, we have 

\T(u)-T(v)\ x 

< \Mu) - Ji(u)|* + |J 2 (n) - J 2 (v)\ x 

^ 2 f\ \2 I |2 ill 2 I I 2 \ * I I (A -\n\ 

S C 2 ^Plc([0,T];if) ' l U lL 2 (0,T);V + l v IC([0,T];tf) ' \ V \L 2 (0,T);Vj ' \ U ~ V \X l 4 ' iZ J 

+ {c 4 T)^\u-v\x (by Lemma [42] and |43]) 

1 

< (c 2 M)* (\u\ 2 L 2 {0)T) . v + \v\ 2 L2{0yT) . v y \u-v\x + (c 4 T)^\u-v\ X - 

Due to the absolute continuity property of Bochner integral, we can choose r £ (0, 1] 

s.t. 

"'\L 2 (0,t);V "t" r y lL 2 (0,r);l/ 

Let T = min{r, 1, ^} and X To := C([0, T ]; F) n L 2 (0, T ;V). We have 



"iWrfcV + \v\h(o,TY,vY < (2Mc 2 )-i (4.13) 



|,F(u)- Hv)\x To <{^ + ^)\u-v\x To =-\u-v\x To - (4.14) 

Applying the modified fixed point lemma I4~4"l equation 

u = (f) + F{u) = 5(-)u + z + Ji(u) + J 2 («) (4.15) 

has a unique solution u in C([0,To]; -ff)nL 2 (0,To; V) and the solution satisfies \u\x T < 
M. D 

Secondly we give a priori estimates and obtain the global existence. Denote by u, 
the local solution of (|4.ip over [0,To]. Let v{t) = u(t) — z(t). Then v(t) is the mild 
solution of equation 

v(t) = S(t)u - [ S(t- s)B(v(s) + z(s))ds - [ S(t- s)N(v(s) + z(s))ds. (4.16) 
Jo Jo 

According to [12J 2.1.20, v(t) is the weak solution of the following differential equa- 
tion with random parameters: 

f v(t) 
4/ + Av{t) + B(v(t) + z(t)) + N(v(t) + z(t)) = 0, 
dt ( 417 ) 

_ v(0) = u 

Inspired by [24J Chapter 15.3, we take advantage of the relation between weak and 
mild solution, give a priori estimate which ensures the global existence of solution. 
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Proposition 4.6. J2J/ Assume that v is the solution of (|4,16|) on the interval. Then 
we have 

n C5 L \z(s)\ 2 ,ds „ f eg f UMI 2 1 dr 

sup K*)| 2 < e J ° M \ Uo \ 2 + / e ° Js ^o 1 gi ( s )ds, (4.18) 

tG[0,T] JO 

/ |u(t)|ydt < c g \uq\ 2 + c 5 c 6 sup \v(t)\ 2 \z(s) \ 2 H i ds + c 6 / gi(s)ds (4.19) 
JO te[o,T] Jo ° Jo 

where C5 and C6 are positive constants depending on X\ and O, g\ is an integrable 

function depending on z. 

Since z £ C([0,T]; V), the following theorem is an immediate consequence of theo- 
rem U3] and proposition 14.61 



Theorem 4.7. For all T > and uq € H , the equation (|2.19|) P-a.s. has a unique 
solution u <G C{[0,T}; H) D L 2 (0,T;V) in the sense of KB . 

5 Random attractor 

First we introduce the fractional Ornstein-Uhlenback process (fractional O-U pro- 
cess) to construct the random dynamical system. 
Consider the linear stochastic evolution equation 

dz(t) = Az(t) + dB H (t), teM. (5.1) 

By the similar argument in [25] chapter 3, we can construct a unique stationary solution 
Z(t) (the so-called fractional O-U process) such that 

Z(t) = Z(9 t uj) = f S(t- r)dB H (r) (5.2) 

J — 00 

r° 

Z(u) = lim A / S{-r)dB H {r) (5.3) 

n-^co J_ n 

Next we prove 

1 f n 
lim - / \\Z(e t «,)\\ 2 H idt < 2cof3 D (4) ■ C(4). (5.4) 

where cq is a constant depending on the space injection. In order to convert this 
integration with respect to time variable into a integration with respect to sample space, 
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we may use ergodic theory. Consider the real-valued continuous function \Z(6.tu)\ 2 m , 
we have 



E|Z(w)|^i=E 



E 



n=l 
oo 



S(- 


-r)dB H {r) 


2 


oo „ 

J2 S(-r)e n dp*(r) 

n=l J ~ t 


E 


fs(-r)e 


n dff(r) 



H' 



nh 



zZ}™^ \ K h( S (— ) e n)\ 2 L*(-t,0;Hi) 



n=l 
oo 



V lim \S(— )e„| 2 i_ H 

* — * £— »oo 7-2 



-11 



n=l 

00 



(^(-tAffo 1 )) 



(Since W = J| (^ 2 (-t, 0; #<}))) 



< N^ lim Co |»S'(— -) e nli2(_ t o-h 1 ) (by Hardy-Littlewood theorem) 
n=l °° 

Elim c / \Aie~ tA e n \ 2 dt 
i— »oo / n 
n=l ,yu 

= Y^ lim c / |An< 

*— ' t-i-oo / n 
n=l JU 



e~ tA "e n rdi 



c V lim A„*(l-e-' A ") 

* — » t— >oo 
n=l 



2 C0 /3 D (4)-C(4) 



(5.5) 



Hence, we have iZffl.w)!?,! £ L 1 (il,P). Since (fi, J 7 , {#(£)}*£«) is the metric dynamical 
system, we can use the Birkhoff-Chintchin Ergodic Theorem to obtain 

1 f n 

lim - / \\Z(9 t u)f Hl dt = E\\Z{u)\\ 2 Hl < 2c (3 D (4) • C(4). (5.6) 

n— >±oo n Jn "0 

According to Theorem 14. 5( Vto G R, u(t,u>;to,Uo) is the unique solution of the 
equation 

u(t;t ) = S{t-t Q )uQ- J S(t-s)B(u(s))ds- I S(t-s)N(u(s))ds+ I S(t-s)dB H {s). 

(5.7) 
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In this section, let u(t,uj;to) = v(t,uj;to) + Z(t,ui), we have 
v 2 (t)+ I S(t-s)dB H (s) 

J — oo 



oo 

t 



=S(t)u - S(t-s)B(v 2 (s) + Z(s))ds- S(t- s)N(v 2 (s) + Z(s))ds + S(t - s)dB H (s). 
J to Jto JtQ 

(5.8) 

Since 

f S(t- s)dB H (s) = S(t - t )Z(e to uj), (5.9) 

v(t,uj;to,UQ — Z(9t Q u})) is the unique solution of the integral equation 

v(t) = S(t)(u -Z(9 to u;))- f S(t-s)B{v(s) + Z(s))ds- f S(t- s)N(v(s) + Z{s))ds. 

Jto Jto 

(5.10) 
According to [12] 2.1.20, v is the weak solution of the following differential equation 

nil 

^ + A(v + Z) + B(v + Z) = (5.11) 

v(t ) = u - Z(6 to uj) (5.12) 

We can now define an continuous mapping by setting 

4>(t,uj,uo) =v(t,uj;0,u - Z(uj)) + Z(6 t u), V(i,w,« ) e M x fi x JT. (5.13) 

The measurability follows from the continuity dependence of solution with respect to 
initial value, the cocycle property follows from the uniqueness of solution. Thus, <f> is a 
RDS associated with f|4.1|) . 

In the rest of this section, we will compute some estimates in spaces H and V. 
Then we use these estimates and compactness of the embedding V ^-> H to obtain 
the existence of a compact random attractor. Assume that C\ is the constant which 
satisfies the inequality of trilinear form b 

b(u,v,w) < CM 1 ' 2 ■ M 1 / 2 • \v\ H i ■ \w\ 1/2 ■ H 1 / 2 (5.14) 



as in 



Lemma 5.1. If cqC 2 < g (4^(4^ > then there exist random radii Ph{u) > and Pi(lo) 
such that for all M > there exists t2{oo) < —1, such that whenever to < t 2 and 
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|uo| < M, we have 

\v(t,uj;to,uo-Z(6 to u;))\ 2 <p H (u;), We [-1,0] (5.15) 

\u(t,u;;t ,uo)\ 2 < p H (u), Vt€[-1,0] (5.16) 

\v(t)\vdt < Pi (oj) (5.17) 



i 
o 

2 



\v{t) + Z{t)\ z v dt < Pi (uj) (5.18) 

l 

Proof. The proof is similar to [8j Proposition 3.7. Multiple (|4.17|) by v(t) and then 
integrate over O. We have 

\'^f- + W)\v = ~Hv(t) + Z(t),v(t) + Z(t),v(t)) - < N(v(t) + Z(t)),v(t) > 
< \b (v(t) + Z(t), Z(t),v(t) + Z(t)) |- < N(Z(t)),v(t) > . 

(5.19) 
The above inequality take advantage of < N(v),v > > (see [3]) and (|2.15p . 
In the sequel we omit the time variable t. Firstly we estimate trilinear form b. 
b(v + Z,Z,v + Z) 



< C\ \V + Z\ ■ \Z\rrl ■ \v + Z\ TTl 

— j- i ii in () i \n 

C\ , l9 ,0 ClC 2l ,n ( 5 - 20 ) 

< -^- Z 2 „i ■ u + Z 2 + -4-^ u + Z 2 Hl 
2C 2 H o 2 ' '#(> 

^^l| 7 |2 |„|2 , /nf /nf |„|2 i ^1| 7 |2| 7 |2 _i_ r> r» I 7l 2 

— 77" Z \N 1 \ V \ + ^1^2P ifl + TT \ Zj \ Z W 1 + W^2 -^ h-1) 

where C 2 is a positive constant which will be specified later. Secondly we estimate 
nonlinear term N. For all n > 0, we have 

-<N(Z),v>< f i e- a / 2 \Z\ H i\v\ H i 

u 2 (5 ' 21) 

<ri\v\%i + -r^-\Z\ 2 Hl . 
M o 4rie a H o 

Comprehensively, 

1 d o "1 i i2 1 i i2 

- — |f| H \v\ H — |v|y 

?T 2 2 r „2 (5-22) 

<^|z|^m 2 + (dc 2 + ri )\v\%r + ^-\z\ 2 \z\ 2 Hl + c x c 2 \z\\, + -^-\z\ 2 Hl , 

where Ai > 4 is the first eigenvalue of operator A. Let g 2 = -^\Z\ 2 \Z\ 2 Hl +C\C 2 \Z\ 2 Hl + 



2 



W 



4rie c 



\Z\ 2 TTl . Then we have 
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By assumption we can choose C 2 G ( coCi/3d(4)£(4), cT ) an d ?"i small enough and we 
have 

By Gronwall inequality, when t G [— 1,0] and to < — 1> we have 
K*)| 2 

<|u(t )| 2 e V / + / <te(si)e V / tfoi (5.25) 



4( 


*i_ 


Cil^Wl^i' 
c 2 


Ids 
' + 


ft "A ( 

/ S2(si)e 

'to 


hL- 


c 2 


\ds2 


£( 


^1. 


OiUMI^ 

c 2 


Ids 
' + 


/"0 - J^ i 

/ 92(si)e 

ho 


(* 


Ci|z( S2 )l^i' 

c 2 


)ds 2 



<|v(t )| 2 e v / + / #2(si)e v / dsi 

•/to 

Due to the ergodic property of fractional O-U process ([5.6p . we have 



lim 

to^-00 -t Jto 

Choose r 2 small enough such that 



1 f° 

— / \Z(s)\ 2 „ids = E\Z(co)\„i. (5.26) 

-to Jto H ° ° 



^E|Z( W )| ff a < ^2c /3 D (4) • C(4) < 2 - r 2 < ^ - r 2 . (5.27) 



Then there exists ti(u>) < — 1, such that when to < ^1 we have 

'to 

By Lemma 2.6 of [25], 52 has at most polynomial growth as to —> —00 for P-a.s. 
lo € $7. Thus, we have 



r 

b(t)| 2 <e (1+io)r2 |u | 2 + / e {l+to)r2 g 2 (s)ds, Vie [-1,0]. (5.28) 

Jto 



'to 
fO 



/■0 />0 

/ 5 2 (s)e (1+s)r2 ds < / 52 (s)e (1+s)r2 (is < 00, P-a.s. 

Jto J —00 



(5.29) 



Let p H = 4,j_ oo g 2 (s)e ( - 1+s ^ 2 ds + 2sup fe[ _ 1)0] \Z(t)\ 2 and there exists t 2 (w) < ti(u) < 
— 1 such that for all \uq\ < M 

\v(-l,co;t ,uo - Z{e to uj))\ 2 < 2 / 52 (s)e (1+s)r2 ds (5.30) 



\u(-l,u>;to,u )\ <2\v(-l,ui;t ,uo- Z(6 to ui))\ +2 sup |Z(t)| 

^I- 1 - ! (5.31) 

<p H (v), Vt <t 2 ,tG [-1,0]. 
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In the following we consider a bound of J_ 1 \v(t)\ydt. Integrating (|5.23p over [—1,0] 
we have 

K0)| 2 - \v(-l)\ 2 + c, 1 J° \v(t)\ 2 v dt < f 92{t)dt + £ (^-\Z(t)\lMt)\ 2 dt. (5.32) 

When to < ^2 we have 

o r o r< „ r® 

2 J4. / „ r / „ U\M , ^ l P H / I7W|2j, , l„./ 1\|2\ a 



\v(t)\ v dt<c 6 ( g 2 (t)dt + ^^ \Z(t)\idt + \v(-l)\ z )^C(uj). (5.33) 

1 J-l ^2 J-l 

Similarly, 

I \v(t) + Z(t)\ 2 v dt<2c 6 (J g 2 (t)dt + ^^J \Z(t)\ 2 dt + 2J \Z(t)\ 2 v dt)±C(Lu). 

(5.34) 
Let pi{oo) = max{C(u;), C(u>)} and the proof is complete. □ 

By the same argument as in [8] Lemma 4.3, we have the following lemma. 

Lemma 5.2. Under the assumption of Lemma \5.1l there exists a random radius pv(w) 
such that for all M > and \uo\ < M, there exists £2^) < —1 such that P-a.s. 

\v(t,uj;t ,uo - Z(9 to uj))\l < pv(u), (5.35) 

\u(t, u; t , u )\l <pv{u), Vt <t 2 ,te[--,Q}. (5.36) 

Lemma [5.21 shows that there exists a bounded random ball in H 1 which absorbs any 
bounded non-random subset of H. Since H 1 is compactly embedded in H, we have 
establish the existence of a compact random absorbing set in H. We now state the final 
theorem and the proof is similar to that of Theorem 5.6 in |10| . 



Theorem 5.3. If cqC 2 < „ (4 w 4 n , then the random dynamical system associated with 
(|4.ip has a random attractor. 
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